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^ Abstract 

The influence of the quadrupole shape fluctuations on the dipole vibrations in transitional nuclei 
00 is investigated in the framework of the Instantaneous Shape Sampling Model, which combines 

the Interacting Boson Model for the slow collective quadrupole motion with the Random Phase 
Approximation for the rapid dipole vibrations. Coupling to the complex background configurations 
is taken into account by folding the results with a Lorentzian with an energy dependent width. The 
low-energy energy portion of the 7- absorption cross section, which is important for photo-nuclear 
processes, is studied for the isotopic series of Kr, Xe, Ba, and Sm. The experimental cross sections 
are well reproduced. The low-energy cross section is determined by the Landau fragmentation of the 
dipole strength and its redistribution caused by the shape fluctuations. Collisional damping only 
wipes out fluctuations of the absorption cross section, generating the smooth energy dependence 
f — , observed in experiment. In the case of semi- magic nuclei, shallow pygmy resonances are found in 

agreement with experiment. 



PACS numbers: 21.60.Fw, 21.60.Jz, 23.20.Lv, 25.20.Dc, 27.60.+j 

ON 

o 
o 

> 



1 



I. INTRODUCTION 



The cross sections of photo-nuclear processes, such as (7, n), (n,j), (j,p), {p, 7), (7, a), 
(a, 7) are key elements in various astrophysical scenarios, like supernovae explosions or 7-ray 
bursts. Precise values of these basic data are also indispensible for simulations of processes 
of nuclear technology. Many of the relevant reactions involve unstable nuclei for which 
measurements of the cross sections are not possible. Therefore, the cross sections have to 
be taken from theory. In many cases the statistical Hauser-Feshbach model is applicable, 
which decomposes the total reaction cross section into a product of the absorption and 
emission probabilities of the particles and 7-quants. Theoretical models that predict the 
dipole strength function for 7-absorption or emission through the whole nuclear chart are 
therefore of utter importance. The reactions take place in an energy interval of a few MeV 
around the particle emission thresholds. Aside from the mentioned applications in nuclear 
astrophysics and nuclear technology, the understanding of the mechanisms that determine 
the structure of the dipole strength function on the low-energy tail of the isovector giant 
dipole resonance is a challenge of its own to nuclear theory. 

The present article proposes and tests a new approach, which we call Instantaneous Shape 
Sampling (ISS). This approach combines the microscopic Quasiparticle Random Phase Ap- 
proximation for dipole excitations with the phenomenological Interacting Boson Approx- 
imation for a dynamical treatment of the nuclear shape. The ISS approach aims at the 
microscopic description of the dipole strength function of the many transitional nuclei rang- 
ing between the regions of spherical and well deformed shapes, which execute large shape 
fluctuations. The special focus of our calculations is the behavior of dipole strength in the 
energy range of few MeV around the particle-emission threshold. These energies are most 
important for the mentioned applications. It is also the range where there are great gaps 
in theoretical and experimental knowledge about the dipole strength function. ISS was sug- 
gested in Ref . p] , where it was applied to the Mo, Zr, and Sr isotopes. Further it was applied 
to 139 La [2J. The present work describes the model in detail and presents additional results 
of systematic calculations in various regions of the nuclear chart. Our paper is organized 
as follows. In section [XT], basic features of the photo-nuclear absorption cross section are 
recalled. The ISS is introduced in section III Section [TV] presents the calculational scheme 
of ISS. The version of the Quasiparticle Random Phase Approximation used in this paper is 
laid out in section |V| Section VI contains the detailed discussion of our results. The range 



of validity of our ISS approach is discussed in VII Conclusions are drawn in section VIII 



II. THE ABSORPTION CROSS SECTION 

The cross section <jei{E) for the absorbtion of electric dipole (El) radiation at the energy 
E by an even-mass nucleus is 

tr m {E)= 4M7ES E1 (E), S E1 (E) = ^ (1) 

where the strength function Sei{E) is the derivative of the reduced transition probability 
B(E1; E) f for a transition from the + ground state to a 1~ excited state at energy E. The 
units in Eq. ([!]) are E in MeV, a E \ in mb and B(E1) in e 2 fm 2 . For magnetic dipole (Ml) 
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*{E,r) = <T R — , m L,L^ r ,2 (3) 



radiation one has 

a Ml {E) = 0.0452 ES Ml {E) , S M1 (E) = B{M ^ E)t , (2) 

where the strength function Smi(E) is the derivative of the reduced transition probability 
B(M1;E) ^ for a transition from the + ground state to a 1 + excited state at energy E. 
The units in Eq. ^ are E in MeV, om\ in nib and B(M1) in [i 2 N . The absorption cross 
section enters directly the total cross section for (7 , particle) reactions. For the (particle,^) 
reactions, the dipole strength function determines the 7-cascade depopulating excited states. 
We present only results for the absorption cross sections, because the strength function can 
be easily obtained from the absorption cross sections by means of Eqs. ([Tj [2]). 

The El part of the absorption cross section dominates in the energy range above 6 MeV, 
where the Ml part amounts to only a few percent of the total cross section (see Sec. fV| ). 
Thus, in reviewing previous work we focus on the electric part. The prominent structure of 
the cross section ctei(E) is the Giant Dipole Resonance (GDR). It can be approximated by 
a Lorentzian [3j H] 

E (TE) 2 
E~ R {E 2 -Elf + {EVy 

and is determined by three parameters, the resonance energy Er, the maximum height ctr, 
and the width T. This expression represents the amount of radiation absorbed by a classical 
damped dipole oscillator, where T/h is the dissipation rate of the vibrational energy which 
determines the width of the resonance. Concerning the low-energy tail, the question arises 
if and in which way T depends on the energy. The authors of Ref. [3] suggest r oc E 2 
based on the theory of Fermi liquids, where E x is the photon energy. Other authors consider 
the width as constant. For axially deformed nuclei, the GDR splits into two peaks [21 E], 
for triaxial deformation into three peaks, where each of the modes can be described by a 
Lorentzian [8]. The different resonance frequencies reflect the different wavelengths along 
the principal axes of the deformed nucleus [I] . The Hauser-Feshbach codes, used to calculate 
reaction rates for various applications, like such as the one in Ref. [9], traditionally employ 
some version of the two-Lorentzian model [10J. The deformation parameters of the assumed 
axial shape are taken either from experimental B(E2) values or from calculations by means 
of the Micro-Macro (MM) Method [HUTS] . 

In Ref. [3] the experimental absorption cross section was analyzed in terms of a model of 
three individual Lorentzians with relative resonance energies related to the nuclear deforma- 
tion according to hydrodynamics (see, e.g. Ref. jl]) and T oc Ejf, where Er is the energy of 
the resonance. Recently, cross sections in the GDR region were combined with novel photon- 
scattering data obtained from experiments at the ELBE accelerator ^BTTTfj . which provided 
for the first time cross sections from the low-energy region up to the GDR region. The new 
data suggest that the smooth Lorentzian extrapolation to energies far below the peak region 
of the GDR provides a rough estimate for the average trend of the cross sections, but is 
not capable of describing resonance phenomena observed in this energy region. Although 
the experimental information about the cross section below the particle emission threshold 
is still sparse there is evidence for pronounced fluctuations and resonance-like structures 
that are above the Lorentzian [131 EE HE]- Furthermore, below 5 MeV the cross section 
changes more and more in to a discrete spectrum of individual 1 states. This energy region 
becomes important in reactions involving neutron-rich nuclei. The dipole strength function 
that determines the 7-cascade in (particle^) reactions may also belong to this energy region 
if the Q value is low. 
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FIG. 1: (Color online) Absorption cross section of 130 Ba calculated in QRPA with the equilibrium 
deformation (3 = 0.171 (cf. Tab. [TT]) . Solid black curve: width T = 3 MeV, thin blue curve: width 
r = 0.1 MeV. 

A more microscopic approach is needed to account for these shortcomings. Moreover, a 
microscopically founded description of the dipole strength function is a longstanding chal- 
lenge of its own to nuclear theory. The traditional approach is to start from the Quasiparticle 
Random Phase Approximation (QRPA) [18], which extends the mean field description of 
the nucleus by taking into account the quasi-bosonic part of the residual interaction, which 
generates a coupling between the elementary two-quasiparticle excitations in the selfconsis- 
tent mean field potential [19]. After the QRPA dispersion relation is solved one obtains a 
discrete series of states that are superpositions of the two-qpasiparicle excitations. These 
states describe how the cross section of the collective dipole vibration is spread over the 
two-quasiparticle excitations. This mechanism is called Landau fragmentation (or Landau 
damping in the case of a dense two-quasiparticle spectrum). It accounts for a part of the 
width of the GDR. 

For example in Fig. [TJ we show the absorption cross section for the axially deformed 
nucleus 130 Ba. The discrete energies of the QRPA solutions are folded with a Breit-Wigner 
distribution with a width of V = O.lMeV (see Eq. Q below). The strong fluctuations of 
the cross section reflect the individual structure of the two-quasiparticle excitations that 
contribute to the GDR. The nucleus is prolate. The two main peaks correspond to the 
vibrations along the long axis (K = 0, low) and along the short axis (\K\ = 1, high). The 
Landau fragmentation generates a cross section a at low energy, which, as discussed in the 
following, is of the order of the observed one. However, there must be strong additional 
couplings of the QRPA "doorway" states to more complicated excitations, which smooth 
out the fluctuations completely. 

In standard QRPA calculations these couplings are taken into account in a phenomenolog- 
ical way by folding the QRPA solutions with a Breit-Wigner or Lorentz function analogously 
to the three solutions of the dipole oscillator just mentioned. We prefer for simplicity the 
Breit-Wigner distribution which is practically equivalent to the Lorentz function used in Eq. 
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<J3J) . The Breit-Wi gner distribution is given by 

r /27r 

Y) = °u {Ev _ E) 2 + T 2 /A > W 

where a v is the cross section of the QRPA solution at the energy E v . The cross section takes 
on the smooth shape seen in experiments, if folded with a large width T = 3 MeV. In the 
considered case the relation Er(K = 1) — Er(K = 0) < F holds, and the two peaks merge 
into a broad peak, the large width of which reflects the deformation. QRPA calculations for 
static spherical or deformed shapes (see, e.g. Ref. [201 - 122] ) use such a smoothing procedure 
to account phenomenologically for the neglected coupling to more complex configurations. 
Accordingly, in Refs. [231 121] strength functions for applications in Hauser-Feshbach codes 
were calculated throughout the nuclear chart. The authors start from the Skyrme Hartree- 
Fock-Bogoliubov (HFB) mean field, which provides the deformation parameters. However, 
they carry out the QRPA calculations for spherical shapes only and include the splitting 
caused by deformation in a phenomenological way. Hence, QRPA describes the parts of the 
GDR width that originate from Landau fragmentation and static deformation. 

The microscopic origin of the remaining part of the spreading of the GDR has been re- 
viewed, for example, in Ref. [25] [26J. The escape width accounts for the emission of 
particles from the QRPA states. It is only of importance for light nuclei or nuclei near the 
drip lines. The spreading width describes the coupling to more complex configurations. 
It is suggested that the coupling to correlated four-quasiparticle states should contribute to 

whereas the coupling to six-quasiparticle and higher-order quasiparticle excitations will 
only wash out some residual fluctuations of the strength function. For spherical nuclei, the 
coupling to the four-quasiparticle configurations has been taken into account in the frame- 
work of the quasiparticle-phonon models, such as QPM (2Ql E3 |28] , QRPA-PC [29J, QTBA 
[30H32] and RQTBA [33] [31]. The resulting strength functions reproduce the spreading of 
the GDR. Although a generalization to nuclei with static deformation is possible [20] , calcu- 
lations have been restricted to spherical nuclei so far because of the substantial increase of 
the numerical work. However, principal problems arise in transitional nuclei that undergo 
large-amplitude shape fluctuations. 



III. THE ISS MODEL 

In the following we suggest an alternative approach. We explicitely describe the coupling 
of the dipole vibration to the two-quasiparticle excitations by means of the QRPA for de- 
formed shapes. Out of the coupling of these QRPA doorway states to the more complex 
configurations we only take into account the low-energy collective quadrupole excitations, 
which represent the softest mode that couples most strongly to the dipole mode. The 
quadrupole mode is described by a model that allows for large amplitude motion, i.e. one 
that is suited for transitional nuclei, which are the main object of our work. The typical 
frequencies hu(2 + ) of collective quadrupole excitations are smaller than 1 MeV, which means 
about a factor of 10 less than the energies ftw{l~) of the dipole excitations. Because the 
quadrupole motion is much slower than the dipole one we use the adiabatic approximation: 
By means of QRPA, we calculate the dipole absorption cross section a EltM1 (E, (3 n , j n ) for a 
set of instantaneous deformation parameters (/3 n ,7n) of the mean field. We determine the 
probability P((3 n , j n ) of each shape being present in the ground state and obtain the total 
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cross section as the incoherent sum of the instantaneous ones, 



a ISS (E) = J2 P (Pn,Jn) [<TEl(E,p n ,J n )+<J M l(E,p n ,J n )} . (5) 
n 

The zero-point motion in the ground state with respect to the collective quadrupole modes is 
represented by the set of instantaneous shapes (/3 n , j n ) and their probabilities P((3 n , j n ). As 
a first step, the "dynamical" ground state is constructed in the framework of the Interacting 
Boson Approximation (IBA) [35]. IBA is not a compulsory choice. Other approaches that 
describe large-amplitude quadrupole motion could be used as well. However, IBA treats the 
large amplitude quadrupole motion in an efficient way and has proven to be successful for 
a systematic description of transitional nuclei. In practice it is very easy to handle because 
it has a minimum set of free parameters. Moreover, it allows us to generate the discrete 
set of instantaneous shapes (/3 n ,j n ) and the probabilities P(/3 n ,j n ) in a simple way. In the 
next step, a series of QRPA calculations is performed, where (/3 n ,7n) defines the shape of 
the Woods-Saxon (WS) potential in the QRPA Hamiltonian. Then the total cross section 
is obtained as the incoherent sum Eq. (|5p of the respective cross sections <jei(E, /3 n , 7„) 
and ctmi{E, n , j n ) multiplied by the probabilities P(fl n ,j n ). Finally, the coupling to the 
more complex configurations is taken into account by folding (Jiss(E) with the Breit-Wigner 
function 

Wd ( £ ) = / ^a ISS (E-) {E _ r ^ (E , m (6) 

The width is chosen to depend quadratically on the photon energy E as expected for colli- 
sional damping (CD), i.e. T(E) = aE 2 (5], E3 |26] . We will refer to this phenomenological 
correction as "collisional damping" (ISS+CD) although it comprises all kinds of couplings 
that are not explicitly taken into account. The value a = 0.0111 used in our calculations cor- 
responds to T(E) = 2.5 MeV at E = 15 MeV. The relation hu(2 + )/hu(l-) < 1 is taken as a 
justification for neglecting possible phase correlations between the different shapes (/3 n ,7n). 
In other words, we assume that the deformation does not change during the excitation of 
the nucleus by the absorbed photons, that is, the photon "sees" the shape of the nucleus 
that absorbs it. Thus, the photon current on a target samples the various instantaneous 
shapes of the nuclei in the ground state. Accordingly we suggest the name Instantaneous 
Shape Sampling (ISS) QRPA for the approach. In the following the different steps of the 
ISS-QRPA approach are explained in more detail, results of calculations for a selection of 
nuclei are presented, and a more sophisticated justification of the ISS procedure is given. 



IV. THE CALCULATION OF THE INSTANTANEOUS SHAPES AND PROBA- 
BILITIES WITHIN IBA 

The family of collective quadrupole states in transitional nuclei is described by the sim- 
plified IBA Hamiltonian suggested in Refs. [36j EZ] 

H IBA = c[{l-Qn d -^Q-Q\, (7) 
Qv = s ] d^ + d\s + X [d ] ® d } 2 „. (8) 
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FIG. 2: (Color online) IBA symmetry triangles for the even- A isotopic chains 78_86 Kr, 124_136 Xe, 

130-138 Ba and 144-154 Sm _ 

The operators s*, s and dl, denote the I = and I = 2 boson operators, respectively, of 
the IBA-1 model [35], and is the number operator of d-bosons. The factor c sets the 
energy scale and has no influence on the structure of the states. The isoscalar quadrupole 
operator appearing in the Hamiltonian is also used to define the electric quadrupole (E2) 
transition operator 

T^(E2) = e b Q^ (9) 

where e b represents the boson effective charge in units e b. In IBA applications the total 
number Nb of s and d bosons is usually chosen to be half of the number of valence particles 
or holes relative to the nearest closed shells in the considered nucleus. In contrast, we fix the 
values to Nb = 10. The reason is that Nb determines the number of instantaneous shapes 
to be sampled. The small values of N b near closed shells would not allow us to represent 
the fluctuations of the shape of these nuclei. Fitting the IBA parameters, we found that 
we could reproduce the low-energy states with the same accuracy as for the traditional IBA 
choice of N b . The range of the essential parameters ( and x in ([7j [8]) is restricted to the 
intervals [0,1] and [0, - a/7/2 ], repectively. Within this range, which defines the so-called 
symmetry triangle shown in Fig. [2j the whole variety of transitional structures between the 
vibrational, rotational and 7-independent nuclei is included. For a given nucleus the values 
of ( and x are searched by a fit to the experimental energy ratios E(4i)/E(2i), E(02)/ E(2i), 
and E(2 2 )/E(2 1 ). The resulting IBA parameters values (C, x) f° r f ne nuclides considered 
in this paper are collected in Tab. [T| Fig. [2] places the values for the Kr, Xe, Ba and Sm 
isotopes into the symmetry triangles, where the definition of the polar coordinates (p, 8) 
according to Ref. is applied. Such contour lines within the symmetry triangle may help 
to extrapolate to nuclei with less well known excitation spectra. 

The next step is the calculation of the probability distribution P{fi n , r ) n ) of the deforma- 
tions P> n and 7 n in the + ground state of the boson Hamiltonian Eq. ([7]). We follow the 
method proposed in Refs. [381 EH]- We relate the IBA to the shape of the Woods-Saxon 
potential used in QRPA by assuming that its deformation parameters /3, 7 are the same as 
in the expression for the electric quadrupole (E2) transition operator of a charged liquid 
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TABLE I: Optimal IBA parameters £, x an d the boson charge for the considered nuclides. 



A 

Kr: 78 
80 
82 
84 
86 



JV 
~42~ 
44 
46 
48 
50 



c 

0.47 
0.35 
0.0 
0.35 
0.0 



_X 

-0.79 

-0.06 

-1.2 

-0.76 

-1.2 



e(, in e b 
0.088 
0.075 
0.067 
0.044 
0.049 



Xe: 124 
126 
128 
130 
132 
134 
136 



70 
72 
74 
76 
78 
80 
82 



0.63 
0.61 
0.63 
0.60 
0.53 
0.0 
0.0 



-0.04 

-0.01 

-0.01 

-0.03 

-0.0 

-1.20 

-1.20 



0.094 
0.086 
0.083 
0.079 
0.074 
0.082 
0.085 



Ba: 130 
132 
134 
136 
138 



74 
76 
78 
80 
82 



0.61 
0.60 
0.55 
0.0 
0.0 



-0.26 
-0.03 
-0.01 
-1.20 
-1.20 



0.105 
0.091 
0.084 
0.091 
0.068 



Sm: 144 
146 
148 
150 
152 
154 



82 
84 
86 
88 
90 
92 



0.0 
0.57 
0.54 
0.64 
0.62 
0.71 



-1.20 
-0.02 
-0.59 
-0.41 
-1.15 
-1.31 



0.072 
0.071 
0.087 
0.108 
0.146 
0.142 



drop, 



3ZeR< 

4:71 



P 



D 2 



cos 7+ (D% + Dl 



sin 7 



(10) 



where the Df„, functions transform the quadrupole moments from the body-fixed frame to 



[IV 



the laboratory frame in the standard way. From Eqs. (|9| and (10) it follows that the two 
scalar invariants constructed from the IBA quadrupole operator given by Eq. (|8~1) are 
connected with the deformation as follows 



Q2 
<Z3 



[Q 
[Q 



Q} oc (3\ 

[Q®Q} 2 }r ) oc/3 3 cos37. 



(11) 
(12) 



A set of localized states \n) is generated by diagonalizing q 2 and q% within the basis of 
Nf, = 10 of s-d-boson states. Because the scalars q 2 and commute, they can be simul- 
taneously diagonalized. As we are interested only in the probability distribution P(/3„,7„) 
of the IBA ground state |0+), the diagonalization is restricted to the set of + basis states 
within the boson space of maximal N b = 10 d-bosons. The eigenvalues g 2 ,n and g 3>n provide 
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the values of the deformation parameters (/3 n ,7 n ) that are assigned to each localized state 
\n) by the relations 



£ = y/E (*£h) **■ cos37 " = \/^!(^ (13) 

where e& is the effective boson charge and R = 1.2A 1 / 3 fm is the nuclear charge radius 
(cf. Eqs. ^ and (10)). The eigenstates |n) of the operators (11, 12) are identified with 



instantaneous" mean field states with the deformation parameters {J3 n ,j n ). As usual, the 
boson charge is adjusted to measured B(E2, 2+ — > 0+) values , which are taken from 
the compilation [ID] . The probabilities P(/3 n ,7 n ) are finally obtained by projecting the 
eigenstates |n) onto the IBA ground |0 + ), i.e. 

Wn,7n) = |(0+|n)| 2 . (14) 

Our procedure assumes that the instantaneous charge density and the instantaneous mean 
field have the same deformation. 

As an example, we discuss the IBA part of the ISS calculations for the even- A chain 
78-86 Kr. In the upper part of Table [i] the IBA parameters ((, x) are given, and in Fig. [3] 
the resulting instantaneous deformations (/3„,7 n ) and the probabilities P n are displayed. 
Note that (/?, 7) are curvilinear coordinates with the volume element dft^d'y sin 37. These 
geometric factors are part of the probabilities P n . This has the consequence that 7 = 30° is 
favored. For example, 80 Kr has \ ~ 0, which means that it tends to a 7-instabliblity between 
the U(5) - 0(6) limits (cf. Fig. [2|. The corresponding distribution in Fig. [3] looks as if the 
nucleus had a stable triaxial deformation, which however only reflects the factor sin 37 in 
the volume element. In a similar fashion, (5 = is suppressed. For example, 86 Kr has ( = 0, 
which means it is a spherical vibrator (U(5) limit). The corresponding distribution in Fig. [3] 
looks as if the nucleus had a stable deformation, which however only reflects the factor fF 
in the volume element. 



V. THE QRPA 

The Hamiltonian H for the QRPA [19J comprises the mean field part h m f and a residual 
dipole-dipole interaction v res , 

H = h ( ^ f (f3 n , ln ) + h$(p n , 7 n) + v res , 

0/?n,7n) = 4s(/?n,7n) + A W (P +W + P«) - A«iV«, (15) 

where the indices r = p,n refer to protons and neutrons, respectively. The spatial part of the 
mean field hyysiPniln) is the triaxial Woods-Saxon (WS) potential (for parameters see Ref. 
[IT]). The deformation parameters (/3 n ,7n) of the WS potential define the instantaneous 
quadrupole-deformed shape about which the nucleus executes the isovector dipole oscilla- 
tions. In order to describe the partial occupancies of the single-particle levels in open-shell 



nuclei a static monopole pair potential is supplemented to the WS part in Eq. (15), which 
is defined by the gap parameter A^ T ^ and the Fermi energy A (t) . Therein, and 
denote the particle number and monopole pairing operators for the protons and neutrons, 
respectively. 



9 




0.3 " 0.3 




FIG. 3: (Color online) Probability distributions for 78 86 Kr. The instantaneous nuclear shapes 
over the f3 — 7 plane were found by means of the IBA with the parameters given in Table [TJ 



The residual interaction consists of two terms, v res = VreP + viel l \ where the index El 
refers to the negative-parity electric excitations and the index Ml to the positive-parity 
magnetic excitations. The El and the Ml dipole excitations of the QRPA are calculated 
separately, and their contributions to transition strengths are added up. Our QRPA calcu- 
lations assume schematic interactions of the dipole-dipole type. The El part is given by the 
electric dipole-dipole interaction 

vifp = \xtA £ %Y + \^{ £ « *t) 2 , (is) 

i=l,A i=l,A 

where here r = ±1 holds for neutrons and protons, respectively. The interaction Vres 
is expressed in terms of the doubly stretched coordinates x" referring to the nuclear self- 
consistency model (12] . The term v rfP is the simplest ansatz for a residual interaction with 
the signature 1~. The inclusion of an octupole-octupole term has been investigated too. 
It is not taken into account here, because it was found that it had practically no effect on 
the El cross section. It is worth mentioning that according to the investigations [22j S3], 
the isovector dipole term is the most important contribution in an expansion of a realistic 
Skyrme-type interaction into separable interaction terms. As suggested in Ref. jH], the 
isoscalar (t = 0) interaction term is used for removing the spurious center-of-mass motion. 
Choosing a value x t=Q = 1000 MeV fm -2 for the isoscalar coupling constant ensures that 
the QRPA states have no spurious contributions. The isovector strength constant K t= \ 
determines the mean position of the GDR. The A-dependence of the isovector strength K t= \ 
is assumed to be given by the selfconsistent strength factor [32] 



MUJ o 

x t= i = -x sc r] = — 77, (17) 

where M is the nucleon mass and ftw = 41A 1 / 3 MeV. For the remaining factor rj an empirical 
value rj pa 3 fm -2 is estimated from the systematics of the GDR peak energy. In our QRPA 
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FIG. 4: (Color online) Phot-absorption cross section of 138 Ba calculated with ISS-QRPA and 
energy-dependent width T(E X ) = 2.b(E x /lh) 2 MeV. Lower panel (b): E x = 6 - 10 MeV. Red solid 
curve: El contribution, red dotted curve: Ml contribution, solid green curve: summed El + Ml 
cross section. Upper panel (a): E x = 10 — 20 MeV. Black curve: QRPA with the equilibrium 
deformation in Tab. [IIJ green curve: ISS-QRPA. 





calculations, the 77-value is kept constant within an isotopic chain and is adjusted to the 
empirical peak position of the GDR of one of the isotopes in each chain. 

The Ml part of the residual interaction v res is the magnetic dipole-dipole interaction 



res 2 / j ' ~3~ ~ 2 

t=0,l t=0,l 



The terms J* • J* and S* • S* are composed of the isoscalar (t = 0) and isovector (t = 1) parts 
of the total angular momentum operator J = L + S and the spin operator S, i.e. J* =0,1 = 
J (p) + (-l)'J (n) and S^ ' 1 = S (p) + (-l)*S (n) . A possible quadrupole-quadrupole interaction 
term turned out to be unimportant for the Ml strength above 4 MeV excitation energy. 



The same magnetic dipole interaction as in Eq. (18) was used in Ref. |45j to describe the Ml 
properties of the Mo isotopes. In these investigations it turned out that the Ml transition 
strength in the interesting energy region above 5 MeV is completely dominated by spin 
vibrations that are generated by the strong repulsive isovector spin-spin term S* =1 • S* =1 in 
Eq. (18). The related strength parameter k* =1 is not well known from literature. Therefore, 



we have used the same value k^ =1 = -1 MeV/H 2 as in our previous study [15]. The isoscalar 
spin-spin term S* =0 -S* =0 is neglected because the isoscalar spin part S <=0 in the Ml transition 
operator is reduced by about a factor of twenty as compared to the corresponding isovector 
spin part. However, a large isoscalar term J* =0 • J* =0 is included by choosing ^j =0 = 1000 
MeV/h 2 in order to eliminate effects of the spurious rotational motion. The isovector term 
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jt=i _ jt=i - g QU £ b ecause only influences the scissor mode in deformed nuclei, the Ml 
strength of which appears below 5 MeV [45]. This energy region is not considered in the 
present work. 

The Ml strength is essentially generated by spin-flip transitions between high-j spin- 
orbit partners. Accordingly, it is expected to produce a summed B (Ml) strength of a few 
fi 2 N distributed in the energy region 7-9 MeV in medium- mass nuclei (see, e.g. [IE]). This 
corresponds to a cross section contribution of a few millibarn. In some cases this value can 
reach up to 20 % of the corresponding El cross section in the low-energy region. As an 
example, we present the QRPA result for the nuclide 138 Ba in Fig. |4j showing the El and Ml 
contributions separately. Because the measurement of the parity is a quite challenging task 
there are hitherto only a few experiments that identify the Ml part of the dipole strength. 
For this reason we do not display the relatively small Ml contribution to the dipole cross 
section separately. 

Carrying out the QRPA begins with the diagonalization of the deformed WS potential 



of the Hamiltonian (15) where we use an oscillator basis with the shells N = — 8. The 



pair field in Eq. (15) is included in BCS approximation, which transforms the creation and 



annihilation operators and c k of the WS levels k to the quasiparticle operators 



a 



a 



u k ct + v k c- k . (19) 



Here, u k and v k are the usual BCS amplitudes, and k labels the time-conjugate WS levels. 



In terms of the quasiparticle operators the mean field part of the Hamiltonian Eq. ( 15 ) takes 
the diagonal form 

h ™f = Yl £k K a k + a t a - k )' ( 20 ) 



where e k = ■sj [e k — A) 2 + A 2 are the quasiparticle energies in the WS potential. The values 
of the pairing gaps A < - p,Tl * ) are derived from the binding energies by using a five-point formula. 
The standard way of performing the QRPA consists of solving the equation of motion 

[H,nt] QRPA = E u nt (2i) 



for the phonon operators f2+ (cf. Eq. (22) below ) via a matrix diagonalization [IB] . The 



suffix QRPA in Eq. (21 ) indicates that only the quasibosonic part of the residual interaction 



is included in the commutator. The set of eigenvalues E v forms the discrete spectrum 



the QRPA amplitudes (pil and ipul- The phonon operators, 



of the vibrational dipole excitations. The eigenvectors of the Hamiltonian in Eq. (|21|) define 

it and 



^ = E [<& a t a v + V©a*"**], (22) 
kk' 

create the vibrational states v as a superposition of two-quasiparticle and two-quasihole 
excitations. The partial cross section a v for a dipole excitation from the QRPA ground 
state |) to a vibrational state at the energy E v is 

a v {E) = f&E v \(n v M^JI 2 6{E - E v \ (23) 

where M^ o!e means the electric (-) or magnetic (+) dipole transition operator. Mea- 
suring the cross section in mb and the energies in MeV, the respective scale factors are 
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/(-) = 4.037/(e 2 fm 2 ) and /« = 0.0452/(/i 2 v ) (cf. Eqs. Q §)). The total cross section 
<t(.E) is given by summing over all the partial cross sections a u . Finally we replace 5{E — E V ) 
by a Breit-Wigner distribution of finite width T, which gives the previous expression Eq. 
Q for the cross section a(E, F). 

To circumvent the direct evaluation of the equation of motion (21 ) which involves typically 
a large matrix diagonalization with a rank of n pa 10 4 — 10 5 we apply the strength function 
method [20J. With this method, the summation in Eq. Q can be written in terms of a 
contour integral which finally is cast in an analytical formula for the function a QRPA (E,F) 
that is explicitly given in Ref . [2D] . The use of the analytical expression for u QRPA (E, F) leads 
to an enormous simplification of the practical performance of the QRPA which is of crucial 
importance for accounting of the variety of shapes (/?, 7) inherent to the ISS calculations. 
The QRPA calculations are carried out with a constant width of F =0.1 MeV. This width 
is small enough to retain all relevant structure of the cross section. The method is also 
efficient for a separable interaction that contains more terms than the leading dipole-dipole 
term considered in our paper. The authors of Refs. demonstrated that such separable 
interactions well approximated the non-separable interactions of the Skyrme type 



VI. RESULTS 

The results of ISS calculations for the isotopic chains 78_86 Kr, 124_134 Xe, 128_134 Ba and 
i44-i54g m are s j lown j n pig S . [5]- [7j Our ISS studies of 92_100 Mo, 88 Sr, and 90 Zr are published 
in [I]. The calculations for 139 La are published in [2]. The figures show the calculations 
without collisional damping (denoted by ISS) with a constant width T = 0.1 MeV and with 
collisional damping (denoted by ISS+CD), where the latter were obtained by folding the 
ISS cross section with the Breit-Wigner distribution in Eq. ^ with an energy- dependent 
width F(E) = 2.5 (E/15) 2 MeV. For comparison we display also the QRPA results for the 
equilibrium deformations, which are denoted by MM, where we use the values in Table [TTJ, 
which were calculated by Moller et al. [Hi [12] in the framework of their sophisticated 
Micro-Macro (MM) model. 

In order to characterize the strength function in a more global way we calculated the 
moments 



( EU dEE k (S E1 (E) + S M1 {E)) (24) 
Jo 



of our strength functions Sei,mi(E), which are related by Eqs. (Tj2) to the cross sections 



shown in Figs. [5] - [Tj The centroid energies E, the widths 7 and the integrated cross 



sections S presented in Table III are calculated by means of Eq. (24) from the ISS-QRPA 
results including the collisional damping width F(E), which are labelled as ISS+CD in the 
figures. To characterize the low-energy region the energy E u = 11 MeV is considered as an 
appropriate upper integration limit. Concerning the entire GDR region the moments are 
calculated using E u = 25 MeV as upper limit. The interval E = 0-25 MeV contains more 
than 85% of the Thomas-Kuhne-Reich (TRK) sum rule. Using MM-QRPA with collisional 
damping we also obtained the same quantities for the static equilibrium deformations, which 



are quoted in Table IV 



Figs. j5|{7]show that the inclusion of CD damping (ISS+CD) eliminates the fluctuations of 
the ISS cross section. However it barely increases the average cross section in the astrophys- 
ically interesting energy range, which was already found in our previous study of the Mo 
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FIG. 5: (Color online) Cross sections for the isotopes 86 Kr. Left panel: E x = 7-12.5 MeV, 
right panel: E x = 12-20 MeV. Red curve (ISS): calculated with ISS-QRPA (T=0.1 MeV). Dotted 
blue curve (MM): QRPA (r = 0.1 MeV) with the equilibrium deformations in Tab. |n| Green 
curve (ISS+CD): ISS-QRPA averaged with energy-dependent width T{E X ) = 2.5{E x /lh) 2 MeV. 
The arrows labeled by "n" mark the position of the neutron-emission threshold of the respective 
isotope. 



isotopes pQ. Thus the dipole strength in this region results from the Landau fragmentation 
and the ISS fragmentation, which accounts for the various deformed shapes in the nuclear 
ground state. 

In all considered chains there is a structural change from spherical to deformed shapes or 
reverse. The purely collective models for the GDR predict that an increasing deformation 
leads to an increasing low-energy dipole cross section, because the GDR splits into two 
separate peaks pfl, EHHJ IH|. In Ref. |Sj the same relation between deformation and low- 
energy dipole strength was found for the Mo isotopes within the MM+QRPA method based 
on a Nilsson potential, which was confirmed by ISS-QRPA calculations based on the WS 



potential pp. This tendency is also clearly seen for the Sm chain (cf. Fig. [7] and Tab. III). 

According to the equilibrium deformations, the lightest Kr isotope is oblate and all heavier 
ones have a nearly spherical shape but are expected to be soft against deformation. The IBA 
parameters show a trend from 7-unstable (x = 0) to spherical (£ = 0), which is also visible in 
the deformation distributions in Fig. [3j Although Fig. [5] indicates a certain narrowing of the 
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TABLE II: Ground-state deformation parameters from Ref. [T2] for the Kr, Xe, Ba, and Sm 



isotopes. 


A 


P 


7 


Kr: 78 


0.232 


60 


80 


0.062 





82 


0.071 





84 


0.062 





86 


0.053 





Ba: 130 


0.171 





132 


0.158 


20 


134 


0.132 


30 


136 


0.0 





138 


0.0 





Xe: 124 


0.208 





126 


0.170 





128 


0.184 


25 


130 


0.158 


30 


132 


0.0 





134 


0.0 





136 


0.0 





Sm: 144 


0.0 





146 


0.0 





148 


0.161 





150 


0.206 





152 


0.243 





154 


0.270 






GDR with increasing A (cf. right panel), this is not reflected by a decrease of the low-energy 
cross section in the left panel and the £n values in Tab. |III[ which are nearly constant. 
In fact, the S n values for the equilibrium deformation in Tab. |IV| increase. There a two 
mechanisms to explain this unexpected behavior. One is the A-dependence of the GDR 
peak energy, which increases as A~ 1//3 along the isotopic chain. The other can be traced 
to the emergence of resonances, which reflect the bunching of particle-hole excitations due 
to the progressive degeneracy of the single particle levels with decreasing deformation. The 
conspicuous example is the strong peak near 10 MeV in the MM calculations, seen in the 
left panel of Fig. [HJ which carries a summed strength of 22 mb MeV for 80 Kr and 26 mb 
MeV for 86 Kr. The shape fluctuations in the ISS calculation broaden it progressively with 
decreasing A, which reflects the increasing probability of deformed shapes. A flat bump 
remains of the resonance if CD is included. 

According to Table [IT], the Xe and Ba isotopes change from prolate through triaxially 
deformed to spherical shape with increasing neutron number, i.e. the deformation decreases 
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with the mass number A not only in the Kr but also in Ba and Xe isotopes. From macroscopic 
approaches that take into account the splitting of the GDR into two or three peaks caused 
by axial or triaxial deformation, one expects that the values of E 11 also decrease [HHSl HP] - 
This trend is not visible. In all cases the values of X 11 increase in the MM calculations. 
According to the IBA these isotopes tend to 7-instability as seen in Fig. [2] ISS also gives 
an increase of S 11 with A. Similar to the Kr isotopes, the MM results for Xe and Ba 
show strong two-quasiparticle peaks for spherical shape, which are washed out by the shape 
fluctuations in the ISS results. Only the Sm chain, which spans the region from spherical 
to prolate well deformed nuclei, shows the expected increase of S 11 with A. The examples 
demonstrate that the incompletely dissolved particle-hole structures can substantially change 
the absorption cross section near the neutron threshold and can generate an A dependence 
of the low-energy dipole strength that is opposite to the one of the macroscopic approaches. 
Nevertheless, all absorption cross section, except the ones for semi-magic nuclei, smoothly 
increase with energy. 

One notices that for the considered isotopic chains there is only little experimental in- 
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FIG. 6: (Color online) Cross sections for the isotopes 124_134 Xe (left panel) and 130_138 Ba (right 
panel) in the low-energy region 7-11 MeV. Notations are as in Fig. [5j The arrows labeled by "n" 
mark the positions of the neutron-emission threshold of the isotopes. The black dots with error 
bars in 138 Ba display the measured (7, n) cross section from Ref.[23]. The blue dots below the 
neutron threshold are (7,7') data from the recent measurement by Tonchev et al. |5U| . 
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FIG. 7: (Color online) Cross sections for the isotopes 144 154 Sm. Left panel: E x = 7-11 MeV, 
right panel: E x = 7-20 MeV. Red curves (ISS): calculated with ISS-QRPA (T=0.1 MeV). Dotted 
blue curves (MM): QRPA (r=0.1 MeV) with the equilibrium deformations in Tab. [TTJ Green line 
(ISS+CD): ISS-QRPA with a energy-dependent width T{E X ) = 2.5(£ x /15) 2 MeV. The black dots 
are the measured (7,?i) cross sections and the arrows labeled by "n" mark the position of the 
neutron-emission threshold of the respective isotope. 

formation. Concerning the 7-absorption cross sections at higher energies there exist (7, n) 
data for the Sm chain and for 138 Ba [23] but not for the other nuclides. 

The only 7-absorption cross sections for the low-energy region above 5 MeV up to the 
neutron emission threshold exist for 138 Ba from recent experiments with mono-energetic 
7-rays at the HIGS facility [50], which are included in Fig. [6] The data for this A^ = 82 
semi-magic nucleus show a bump at 8 MeV and possibly another at 9.5 MeV. The A^ = 82 
neighbor 139 La was studied in [2]. The data show a bump at 7 MeV and a shoulder around 
9 MeV, which probably represent the same structures as in 138 Ba. The ISS+CD cross 
section for 138 Ba has a shallow shoulder in this region. In case of 139 La [2], the ISS+CD 
cross section has two very broad peaks at 8.5 and 11.5 MeV. Comparing in Fig. [6] the ISS 
curve with the MM curve (zero deformation), one concludes that the structures originate 
from spherical two-quasiparticle excitations, which are strongly fragmented due to the shape 
fluctuations. The deviation of their energies from experiment may indicate that our choice 
of the Wood-Saxon potential does not quite correctly reproduce the single particle levels. 
A peak in the absorption cross section at 9 MeV was also was found in the A^ = 50 semi- 



17 



TABLE III: Integral properties derived from the moments of the strength function Sei{E) 
in Eq. ([I]) as calculated in ISS-QRPA with energy dependent width for the even-A series 
78-86 Kr 130-l36 Ba l24-l36 Xe and 144-I54g m . _ m ^/ m ^ y = ^(m 2 /m ) - (mi/m ) 2 , and the 

integrated cross section £ = 167T 3 / '{2>hc)m\. The respective integration limits E = 11 MeV and 
E = 25 MeV are indicated as an upper index. The percentage of £ with respect to the TRK sum 
rule is given in the columns (%). 



A 


7~p1 1 /TV IT 1 r\ 

h (MeV) 


— 11 n ~\t\ 

7 (MeV) 


1 1 /TV It ~\ T 1 \ 

h (MeV mb) 


(%) 


7~n2f) /TV T t r\ 

h, (MeV) 


7 (MeV) 


L (MeV mb) 


/(V \ 

(%) 


Kr: 78 


9.6 


1.3 


40 


3.5 


16.3 


3.7 


991 


85 


80 


9.6 


1.3 


35 


3.0 


16.4 


3.5 


1024 


86 


82 


9.6 


1.2 


35 


2.9 


16.4 


3.4 


1054 


87 


84 


9.6 


1.2 


39 


3.1 


16.2 


3.3 


1095 


89 


86 


9.5 


1.2 


45 


3.6 


15.9 


3.3 


1141 


92 


Xe: 124 


9.5 


1.3 


97 


5.3 


15.8 


3.0 


1599 


87 


126 


9.5 


1.3 


98 


5.3 


15.1 


3.3 


1725 


93 


128 


9.5 


1.3 


105 


5.6 


15.0 


3.3 


1750 


93 


130 


9.5 


1.3 


109 


5.7 


15.0 


3.3 


1775 


94 


132 


9.5 


1.3 


113 


5.9 


14.9 


3.3 


1806 


94 


134 


9.5 


1.3 


120 


6.2 


14.8 


3.2 


1840 


95 


Ba: 130 


9.4 


1.3 


83 


4.3 


15.6 


3.4 


1746 


91 


132 


9.4 


1.3 


82 


4.2 


15.6 


3.3 


1766 


91 


134 


9.4 


1.3 


85 


4.3 


15.4 


3.3 


1803 


92 


136 


9.4 


1.3 


88 


4.5 


15.4 


3.2 


1835 


93 


138 


9.4 


1.3 


93 


4.6 


15.3 


3.2 


1866 


93 


Sm: 144 


9.4 


1.3 


95 


4.5 


15.8 


2.9 


1886 


89 


146 


9.4 


1.3 


103 


4.8 


15.7 


2.9 


1907 


89 


148 


9.4 


1.3 


110 


5.1 


15.7 


2.9 


1921 


89 


150 


9.5 


1.3 


123 


5.6 


15.7 


3.0 


1923 


88 


152 


9.5 


1.4 


182 


8.3 


15.7 


3.1 


1896 


86 


154 


9.5 


1.4 


223 


10.0 


15.7 


3.1 


1892 


85 



magic nuclides 88 Sr, 89 Y, and 90 Zr [13] - [16] . The ISS calculations [lj give a peak at the 
correct energy, which is also a fragmented spherical two-quasiparticle state. The inclusion of 
CD broadens the structure somewhat too strongly as compared with the more pronounced 
"pygmy resonance" seen in experiment. Also in the case of the = 82 nuclides 138 Ba and 
139 La, the combination ISS+CD seems to damp the spherical QRPA poles somewhat too 
strongly. Although the number of studied nuclei is still too small for definite conclusions, 
one may speculate that the CD width depends stronger on energy than assumed. The 
observation that in many strongly deformed axial nuclei the width of the upper peak of 
the GDR is twice as big as the width of the lower peak might be taken as evidence for a 
strong increase CD with energy. (However, as discussed below, it may be caused by Landau 
fragmentation as well). The quadratic energy dependence adopted in this paper is derived 
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TABLE IV: Same integral properties as in Tab III but here calculated in QRPA with the equilib- 
rium deformations listed in Tab [TTJ 



A 


E u (MeV) 

i-J I 1V1C V J 


^(MeV") 

J I J.VJ.C V 1 


S 11 ('MeV mb) 

£. ' 1 IVIvj V XXX I 


(%) 


E 25 (MeV) 

-L-J I J.VJ.C V J 


J I 1V1C V 1 


S 25 fMpVmb) 

' ' 1 1V1C V XXX U ) 


(%) 


TO- 78 


Q 6 


1 3 

X . (J 


9q 


9 5 


16 5 

1U . O 


3 4 


1000 


86 


80 


Q 7 

•J . I 


1 1 

J. . X 


33 


9 8 

_■ . O 


1 6 6 

XU . u 


3 2 


1 034 


87 

O ( 


oz, 


q 7 


1 9 




9 q 


16 5 

X w . o 


3 9 


i 05q 


87 

O ( 


84 


q 6 


i i 


40 


3 3 


1 6 9 


3 9 


1 008 


80 


86 


9.5 


1.1 


47 


3.8 


15.9 


3.2 


1144 


91 


Xe: 124 


q 5 


1 3 


Qfi 


5 3 


15.2 


3 3 


i6qq 




126 


q 5 


1 3 


Q6 


5 2 


15.2 


3 3 


1 728 


93 

C7 ij 


1 98 


q 5 


1 3 

X . (J 


1 03 


5 6 


1 5 1 


3 3 


1 753 


04 


1 30 

X OU 


q 5 


1 3 

X . O 


1 05 

x \jo 


5 6 

O . U 


1 5 


3 3 


1 778 

X 1 | O 


04 


132 


q.5 


1.2 


105 


5.5 


15.0 


3.2 


1807 


94 


134 


q.5 


1.2 


115 


5.9 


14.8 


3.1 


1841 


95 


Ba- 1 30 


q 4 


1 3 

X . O 


76 


4 


1 5 6 


3 3 


1 750 


09 


1 39 


q 4 


1 3 

X . O 


78 


4 


1 5 6 

X o . u 


3 3 


1 771 


09 


1 34 


q 4 


1 3 

X . O 


81 

O X 


4 9 


1 5 5 

X (J . o 


3 9 

O . Zi 


1 806 


Q9 


1 36 

X *JU 


q 4 


1 3 

X . •) 


83 


4.2 


1 5 5 


3 1 

') . X 


1 835 


03 


138 


9.4 


1.3 


90 


4.5 


15.3 


3.1 


1866 


93 


Sm: 144 


9.4 


1.3 


93 


4.4 


15.8 


2.8 


1890 


89 


146 


9.4 


1.3 


98 


4.6 


15.8 


2.8 


1913 


89 


148 


9.4 


1.4 


108 


5.0 


15.7 


2.9 


1923 


89 


150 


9.4 


1.4 


118 


5.4 


15.7 


3.0 


1929 


88 


152 


9.4 


1.4 


134 


6.1 


15.6 


3.0 


1942 


88 


154 


9.4 


1.5 


151 


6.8 


15.5 


3.0 


1957 


88 



from the schematic model of a Fermi gas. A function with a steeper energy dependence would 
give less CD in the threshold region if its scale is adjusted to reproduce the peak height of 
the GDR. Such a reduction of CD would barely reduce the average absorption cross section 
at these energies (compare the ISS and ISS+CD curves) but give more pronounced pygmy 
resonances. 

In the right panels of Fig. [7j the distance between the two GDR peaks is somewhat over 
estimated for the well deformed isotopes 152 ' 154 Sm. The discrepancy can be traced to the 
large values of the deformation parameter of 0.306 and 0.341, respectively, which reflect 
the large experimental B(E2,2 + — > + ) values [40J. The estimate of the hydrodynamic 
model 2 (E(K = 1) - E(K = 0))/(E(K = 1) + E(K = 0)) = 0.94/3 [U EHSl HZ] gives simi- 
lar splittings of 4.1 and 4.5 MeV, respectively. This indicates some inconsistency between 
the experimental B(E2,2 + — > + ) value and the observed splitting of the GDR. Another 
problem is visible for the well deformed isotope 154 Sm. In experiment, the second GDR peak 
has a somewhat smaller height than the first one, whereas in the calculation it is opposite. 
Since the high-energy peak is two-fold degenerated (K = ±1) it carries twice the strength of 
the non-degenerate (K = 0) low-energy peak. In order to be lower, its width must be more 
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than two times the width of the first peak. However, collisional damping r oc E% gives only 
(16/12) 2 = 1.8. The stronger Landau fragmentation of the low-energy peak acerbates the 
discrepancy in the considered case of 154 Sm. Many of the well deformed nuclei behave in 
the same way: The high-energy peak has the same height as the low-energy one, indicating 
that its width must be about twice [3] . Different versions of the mean field produce different 
pattern of the Landau fragmentation. In the calculations of Ref. [22] for 154 Sm, the Skyrme 
density functional SLy6 produces stronger Landau fragmentation for the K = ±1 peak than 
for the K — peak, resulting in a better agreement with experiment, whereas the function- 
als SkM*, Skl3, and SkT6 give a too high second maximum. For the same nucleus, Ref. 
[19] obtains sufficient Landau fragmentation in the upper peak for the Skyrme functional 
SKM*, however not for SLy4 and SkP. The results for 154 Sm obtained in [22] and [19] with 
SkM* disagree. At this point it seems unclear if Landau fragmentation can account for the 
widths of the two GDR peaks observed in well deformed axial nuclei. 



VII. RANGE OF VALIDITY OF ISS 

The coupling between the low-frequency quadrupole mode and the purely collective GDR 
mode was studied in Ref. [H] by means of the Dynamic Nuclear Collective Model (DNCM). 
The authors find that the dipole strength becomes distributed over several quadrupole exci- 
tations. The collisional damping washes out the discrete structures to a smooth envelop. In 
Ref. [51] the validity of the ISS approximation was investigated in the same model. The dis- 
crete spectrum of the DNCM was compared with the continuous strength function obtained 
by integrating the instantaneous excitation probabilities of the GDR over the probability 
distribution of the shape parameters in the ground state, which corresponds to a dense set 
of sampling points in ISS. The resulting smooth strength function becomes a good approxi- 
mation of the envelope of the discrete lines if 

d[3 Co> 2 

where /3o is the zero point amplitude of the quadrupole vibration and uj\ and U2 are the 
frequencies of the dipole and quadrupole vibrations, respectively. In our preceding paper 
[T], we gave some qualitative estimate that the ISS is applicable to the dipole excitations 
around the particle emission thresholds according to this condition. 

In order to judge the quality of the ISS approximation in this energy region in more 
detail, we studied the coupling of a single QRPA 1~ pole with the collective quadrupole 
(2 + ) mode in a schematic model. Since the 1~ poles of interest are located substantially 
away from the peak of the GDR, it suffices for the following discussion to consider them as 
dressed two-quasiparticle excitations. Accordingly, their energy is not very different from 
the two-quasiparticle energy, and their transition strength is given by the two-quasiparticle 
strength times an effective charge accounting for the screening. (A detailed discussion of 
this approximation is given in |52j.) Thus, it is sufficient to use the energies and transition 
strengths of the two-quasiparticle excitations in the threshold region as a starting point for 
the following schematic model. 

We generated all two-quasiparticle 1~ excitations in the energy interval between 7 and 
8 MeV and studied their dependence on deformation. We followed the development in a 
diabatic way, such that we calculated the overlap of the wave functions of the quasiparticle 
states at adjacent deformations (3 n and (3 n+ i and associated 1~ states by requiring maximal 
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overlap between n and n + 1 (cf. [53]). A typical change of the two-quasiparticle energy 
was found to be 1-3 MeV over the interval 0.1 < fi < 0.3. One such two-quasiparticle state 
is selected and its coupling to an axial collective quadrupole degree of freedom is consid- 
ered, assuming that the coupling is caused by the deformation dependence two-quasiparticle 
energy (as in the DNCM). The quadrupole mode is described by 10 equidistant points 
fi n = 0.04ft, n — 1, .., 10 and 7 = 0. The step size is chosen to roughly agree with the dis- 
tance between the samples in the realistic ISS calculations (cf. Fig. [3]). The Hamiltonian 
for the collective quadrupole motion is given by the matrix 



=h(5 



m+1 



n,m— 1, 



+ V(m)6 ntr . 



(26) 



The first term is the discretized kinetic energy where h determines the inertial mass param- 
eter. The second is the discretized potential. We studied the two potentials 



Vno(n) 
V sw (n) 



D(n-7) 2 , 

for n 
oo for n 



3, 8 
1, 2, 9, 10 



(27) 
(28) 



which are the discrete versions of a harmonic oscillator and a square-well potential, respec- 
tively. The two-quasiparticle energy is taken as E^ 2qp \n) = (7 + Ae (n - 7)) MeV with 
Ae =0.2 and 0.6, corresponding to a change of 1 and 3 MeV over the considered deforma- 
tion range, respectively. The Hamiltonian describing the two-quasiparticle state coupled to 
the quadrupole mode is 



H n ^ = H n % + (7 + Ae(n-7))5 n 



(29) 



Both and H^ 2q ^ are diagonalized numerically. The resulting eigenvalues are E^ c '{i) 
and E^ c ' 2qp \j). The resulting eigenvectors are U n c J and U n c ^ qp \ respectively, where i and 
j label the respective eigenstates. The transition strength from the ground state 1 to the 
mixed excited state j at the energy E^ c,2qp \j) is 



Sex{j) 



(30) 



where s(n) is the 2qp transition matrix element for deformation point fi n . The analogous 
ISS strength for the sampling point n is 



SissM 



s(n)Ui c Mi 



(31) 



which is associated with the energy E^ 2qp \n). Strength functions are generated by folding 
with a Breit-Wigner distribution in order to account for CD, 



S ex {E) = Y J S ex U) 



2vr ((£( c ' 2 «p)(j) - E) 2 + (r/2) 2 ) : 



Siss(E) — .S/ggl 



2ir{{E( 2qp )(n) - £) 2 + (r/2) 2 ) ' 



(32) 
(33) 
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FIG. 8: Ground state eigenvectors for the harmonic oscillator (left) and square well (right). 




2 4 6 8 10 
E/MeV 



2 4 6 8 10 
E/MeV 




2 4 6 8 
E/MeV 



2 4 6 8 10 
E/MeV 



FIG. 9: (Color online) Strength functions of a RPA pole coupled to the collective quadrupole mode. 
Red circles - exact, red curve dashed - exact+CD (r = 0.8 MeV), green squares - ISS, green curve 
- ISS+CD ( T = 0.8 MeV). Change of the pole energy between /3=0.1 and 0.3: upper panels - 1 
MeV, lower panels - 3 MeV. Potentials: left panels - harmonic oscillator, right panels - square well. 



In the following we discuss only the results for s(n) = 1, because the study of non- 
constant s(n) values, which were derived from the two-quasiparticle excitations in the same 
way as the two-quasiparticle energies, led to the same conclusions. The choice h = 1.12MeV 
and D = 0.028MeV gives an excitation energy of E^ c \2) - E^(l) = 0.31 MeV for the 
collective quadrupole mode. The corresponding ground-state eigenvectors, shown in Fig. |8j 
are distributed over 5 sampling points. Fig. |9| compares the exact transition strength (30) 



22 



with the ISS approximation (31). In the case of the harmonic oscillator potential, the 
oscillator length is 3A/3 = 0.12. In the upper panel the two-quasiparticle energy changes by 
1 MeV between /3 = 0.1 and 0.3 (n = 3 and 8, respectively), which corresponds to £ = 2 
(cf. Eq. (25)). As seen, the energies of the coupled states do not agree with the energies 
of the sampling points. Nevertheless, the locations and the widths of the distributions are 
roughly the same. In the lower panels the two-quasiparticle energy changes by 3 MeV 
between /3 = 0.1 and 0.3, which corresponds to £ = 6. Now ISS becomes a rather good 
approximation to the exact transition strength. This agrees with Ref. [51]. Using the 
continuous version of ISS, the authors found that for £ > 5 the strong coupling limit is 
approached and substantiates the discussion in our previous paper pQ. In the case of the 
square-well potential, the strong coupling limit is somewhat slower approached. Presumably 
this reflects the fact that eigenvectors have a lower probability at the turning than for the 
harmonic oscillator. At the turning point, the energy equals the potential energy, which is 
assumed in ISS. 

The strength functions take CD into account by folding the transition strengths with 
a Breit-Wigner distribution. The width T = 0.8 MeV corresponds to the energy-depend 
width at 7 MeV, used throughout this paper. In the upper panels of Fig. |VII[ the energy 
difference between the sample points is much smaller than the damping width, AE/T = 0.2. 
All structure is averaged out. ISS and the exact strength function are practically identical. In 
the lower panels of Fig. VII the results for AE/T = 0.6 are displayed. Some of the structure 
survives. In the case of the harmonic oscillator the ISS and exact strength functions nearly 
agree, because the discrete transition strengths are very similar. In the case of the square- well 
potential the surviving structures disagree. Since the coupling of the two-quasiparticle state 
to the collective mode is weaker, the energies and transitions strengths of the mixed states 
disagree with those of the sampling points, which is transfered to the strength functions. 
Hence the fluctuations of the ISS strength function do not represent physical structures. 
They are just "sampling noise" that should be disregarded. The ISS strength functions 
for most of the nuclides studied in this paper do barely show structure around the particle 
emission thresholds, which means that ISS is reliable. The shallow peaks of the ISS+CD 
strength functions in the Kr isotopes and Sr and Zr isotopes [lj should be considered as 
real structures ("pygmy resonances") predicted by our model, because their width is much 
larger than the CD damping width. 

The number of sampling points is determined by the number of bosons used for diago- 
nalizing the IBA Hamiltonian. In order to study the effects of this coarse graining of the 
collective mode, we decreased the deformation step in our schematic model to A/3 = 0.02. 
The results for the exact strength functions changed only marginally, which means it is suf- 
ficient to restrict the number of bosons to 10. The ISS sampling noise is suppressed on this 
finer grid. All ISS strength function become smooth peaks because AE/T = 0.3. The fine 
structure of the ISS strength function in the left lower panel of Fig. [9] is similar to the exact 
one, because accidentally we chose the sample points near the deformation points where the 
mixed states localize. In the case of a spherical nucleus and/or at sufficiently low excitation 
energy, where CD is weak, one expects observing the mixed states as resolved lines. Ob- 
viously, ISS will not describe these lines individually. Still, the location and the width of 
the distribution of lines (representing sampling points) will correlate with the location and 
width of a fragmented QRPA pole. 
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VIII. CONCLUSIONS 



The method of Instantaneous Shape Sampling (ISS), suggested in our preceding Rapid 
Communication [TJ has been presented in extended form. It relies on the assumption that 
the photo excitation is a fast process as compared to the shape fluctuations of nuclei, such 
that the total 7-absorption cross section is the sum of the absorption cross sections of a 
set of instantaneous shapes, each weighted with its probability being present in the ground 
state. That is, the 7-quant "takes a snapshot of the instantaneous shape of the nucleus" when 
being absorbed. In the present implementation of the ISS concept, the quadrupole motion is 
described by the Interacting Boson Approximation and the 7-absorption is described by the 
Quasiparticle Random Phase Approximation(QRPA) for a deformed Woods Saxon potential 
combined with a dipole - dipole interaction for the El modes and a spin-spin interaction for 
the Ml modes. The ISS concept can be applied to other versions of QRPA and a different 
description of the collective quadrupole mode. 

Studying the coupling between a dipole QRPA solution near the particle threshold and 
the low-lying quadrupole mode in a schematic model, we found that ISS provides a good 
description of the location and the width of the resulting group of levels, although there is no 
one-to-one correspondence between the lines. Taking into account the collisional damping 
(CD) by folding the discrete QRPA solutions with a Breit-Wigner function, the resulting 
ISS+CD strength functions reproduce the exact ones very well, if the distance between the 
energies of the coupled states is smaller than the damping width. 

We applied our version of ISS-QRPA to the chains of the Kr, Ba, Xe, and Sm isotopes, 
which all span the transitional regions between deformed and spherical shape. As in our 
previous study of the Mo isotopes [!}, we find that the dipole absorption cross section in 
the energy region of photo-nuclear reactions is determined by the Landau fragmentation 
and the dynamical deformation. In order to reproduce the broad peak of the Giant Dipole 
Resonance (GDR) additional CD must be introduced, which we assumed to be proportional 
to the square of the photon energy. Its scale turned out to be nearly independent of the 
nuclear mass. CD smoothes out most of the fluctuations of the ISS-QRPA absorption cross 
section, whereby it does not increase the cross section in the energy region of photo-nuclear 
reactions in any substantial way. For all but semi-magic nuclei, the resulting absorption 
cross section increases with energy in a smooth way, as observed. 

Collective hydrodynamic descriptions of the GDR give an increase of the low-energy 
dipole absorption cross section with nuclear deformation. It is caused by the splitting of 
the GDR into a low frequency oscillation along the long and two high frequency oscillations 
along the short axes. In the case of the Mo and Sm isotopic chains, the deformation increases 
with neutron number. ISS-QRPA reproduces the expected increase of the low-energy dipole 
absorption cross section. However, in case of the Kr, Ba, and Xe chains, for which the 
deformation decreases with neutron number, the expected decrease of the low-energy cross 
section is not found. ISS-QRPA predicts a nearly constant value of the cumulative low- 
energy cross section (E x < 11 MeV) when approaching the shell closure. The reason is 
the A -1 / 3 decrease of the GDR peak energy as well as the progressive bunching of the 
two-qpasiparticle excitations when approaching spherical shape. 

In the case of semi-magic nuclei, relicts of these bunches survive the damping by shape 
fluctuations and collisional damping. They appear as broad bumps in ISS-QRPA cross 
section, which may substantially enhance the absorption cross section around the parti- 
cle thresholds. These "pygmy resonances" are two-quasiparticle excitations dressed with 
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isovector dipole vibrations and fragmented by coupling to shape fluctuations. In the case 
of the N = 50 isotones, the position of the resonance is well reproduced. In the case of the 
N = 82 isotones, some discrepancy between the calculated and observed location may point 
to inaccuracies of the single particle levels of the adopted Woods-Saxon potential. 
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